Heavy quark currents in Ultra-High Energy Neutrino 

Interactions 



R. Fiore^t and V.R. Zoller^t 

o 

(N 

Q ■ Dipartimento di Fisica, Universitd della Calabria 

^ • and 

W . Istituto Nazionale di Fisica Nudeare, Gruppo coUegato di Cosenza, 

C^^ \ 1-87036 Rende, Cosenza, Italy 

^ITEP, Moscow 117218, Russia 

^: 

A . Abstract 

^ : 

£j' We discuss heavy quark contributions to the neutrino-nucleon total cross section at very 

high energies, well above the real top production threshold. The top-bottom weak current is 
found to generate strong left-right asymmetry of neutrino-nucleon interactions. We separate 



(N 
> 



\0 [ contributions of different helicity states and make use of the K-factorization to derive simple 



in 



and practically useful formulas for the left-handed (Fl) and right-handed (-Fr) components 
T^ ■ of the conventional structure function 2x^3 = F^ — Fji in terms of the integrated gluon 

t^^ ■ density. We show that Fi ^ Fr and, consequently, xF^ ^ Ft, where Ft is the transverse 

^ . structure function. The conventional structure function F2 = Fs + Ft at Q^ <ti mf appears 

to be dominated by its scalar (also known as longitudinal) component Fs and the hierarchy 
F5 ^ Fi ^ Fr arises naturally. We evaluate the total neutrino-nucleon cross section at 
L^ ultra-high energies within the color dipole BFKL formalism. 
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The neutrino astrophysics [1] deals with neutrinos of uhra-high energy (UHE), well above 
the threshold of the top quark production. The role of heavy quark currents in UHE neutrino 
interactions has been explored extensively [2]. Different approaches to including the top quark 
effects in the UHE neutrino- nucleon cross section cr'^'^ based on [2] are in use [3]. The goal 
of this communication is to show that a specific choice of relevant degrees of freedom greatly 
simplifies the analysis of a'^'^ in terms of the nucleon structure functions (SF). Making use of 
the K-factorization we isolate the leading contributions involving big log's of Q^, mf and 1/x 
to derive simple and practically useful formulas for the UHE neutrino-nucleon SFs. 

The UHE interactions correspond to neutrino energy above Ei^ ~ 10^ GeV. The overall 
hardness scale of the process 

uN -> i^X (1) 

is determined by the gauge boson mass 

g2 ^ ml.. (2) 

Therefore, the W- boson exchange probes the gluon density in the target nucleon at very 
small values of Bjorken x. Consequences of this observation for a'^'^{Eu) are widely discussed 
[4] (for most recent publications see [5]). 

At small X, to the Leading Log(l/a;) approximation, it is legitimate to consider the W- 
nucleon scattering in the laboratory frame in terms of interactions with the target of the qq/- 
pair (color dipole) which the light-cone W-boson transforms into at large upstream distances. 
The dynamics of the log(l/a;)-evolution [6] of the dipole-nucleon cross section is described by 
the infrared regulated BFKL equation with running coupling derived in [7, 8]. 

The differential cross section for the uN interactions is expressible in terms of the scalar 
(also known as longitudinal), left-handed and right-handed structure functions denoted by 
Fs, Fl and Fr, respectively. 



da^^ G^ ( ml, V 

X- — \ vy \ 



2{l-y)Fs + FL + {l-yYFR . (3) 



dxdQ'^ 47r xml. + Q"^ 
The structure functions Fx are related to the absorption cross sections for the W^-boson in 
the polarization state X = S,L, R, 

F,{x,Q') = -^a,ix,Q'), (4) 
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where aw = fl'^/47r, g"^ = G^^^v^/v^ and the W -^ tb transition vertex is 

ght^igv - 9Al5)b. 

In terms of F^, F^ the transverse structure function Ft = F2 — Fs reads F^ = {F^ + -Fr)/2, 
and left-right antisymmetric structure function F3 is as follows 2x^3 = Fl — Fr. Switching 
L to i? in Eq.(3) yields the vN cross section. 

Then, making use of the color dipole factorization [9] we arrive at the leading-log(l/a;) 
K-factorization representation for the absorption cross section for the W^-boson in the polar- 
ization state A 

axix.Q'^) = f dz f £r\^xiz,r)\'^a{x,r) 

1 ^2 

1^ dzld'kl-^asiq'mx,,K'){S,^l + P,^l). (5) 



aw 



We address first the second line of Eq.(5) ^ to derive some useful analytical expressions 
for F\{x, Q'^). The first line of Eq.(5) will be used in our numerical studies of (j^^{Ey). 

No restrictions on the quark transverse momentum, k, and the gluon transverse momen- 
tum, K are imposed in Eq.(5), where 

2n _ dG{xg,K, ) 

O log K,^ 

and G{xg, k^) = Xgg{xg, k^) is the integrated gluon structure function. It is taken at 

, Q' + "''' (6) 

and the transverse mass of the tb pair is 

^2 _ "^' + k' ^ ml + (k - k)^ 
z 1 — 2; ' 

where z is the fraction of the W^-boson light cone momentum carried by the t-quark. 



^To compare Eq.(5) with Eq.(2) of Rcf. [10] one should make the substitution: 
The numerical factor 16 in Eq.(2) of Rcf. [10] should be understood as 4. 



The terms proportional to Sx and Pa describe interaction of the quark-antiquark states 
with the angular momentum L = (S-wave) and L = 1 (P-wave), respectively. The S-wave 
and the P-wave factors in Eq. (5) for the right-handed and left-handed H^-bosons are as follows 

[11] 

Sr = {gv [(1 - z)mt + zmk] + Qa [(1 - z)mt - zmb]f 

Pr = {gv - 9Afz^ + {gv + gAf{i - zf 

Sl = {gv [(1 - z)mt + zrub] - gA [(1 - z)mt - zmh]f 

Pl = {gv - gAf{l - zf + {gy + gAfz" (8) 

and for the scalar/longitudinal polarization [12, 11] 



Ss = {gv/Q') [2Q^z{l -z) + {mt - m,) [(1 - z)mt - zm,] 
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HgVQ^) [^Q^z{l -z) + {mt + nih) [(1 - z)mt + zmf,]f 

Ps = {gv/Q'){mt - m,f + {gl/Q^){mt + m.f (9) 

In the charged current interactions gA = gv = ^- In the neutral current neutrino interactions 
fT^q = iTT-q' and correspoudlug vector and axial-vector couplings are given by the Standard 
Model. 



In Eq.(5) 



$0= ^ '- 

/k k — k\ 



and 



e^ = z{l - z)Q'^ + (1 - z)m^t + zml (11) 

Then, making use of the technique developed in [13] for electro-production processes, separate 
the K^-integration in (5) into the soft, 

K,^ <k^ = e^ + k'' (12) 

and hard, 

k'>W, (13) 



regions of the gluon momentum. For soft gluons, in the DGLAP [14] region, upon the az- 
imuthal integration we get 



^V'^2 ., 2^^ + (k- 



J 2n ^ 



2\2 



(k2+£ 



2^4 



dip^o 2K2k2 ,^^, 

_!_$^ ~ (14) 

2% ° (k2+£2)4 ^ ) 

In Eq.(5) the QCD running couphng as{q'^) enters the integrand at the largest relevant 
virtuality, g^ = max{A;2, k^}. For soft gluons one can take q^ = k'^ and we arrive at the fully 
differential distribution of the t-quark in z and k, 



dax{x,Q 



2^ 



r-^ 



awOisik^)G{xg, P) 



^ ^(k2 + £2)4^ ^(k2+£2)4 

where k'^ comes from Eq.(12). Integrating over k yields 



(15) 



dax{x,Q'^) 2iTaw ( S\ 



i. - 3.^ y + A)".(aG(....=). (16) 

The leading contribution to Eq.(16) comes from k2 < e^ j?). We factored out the slowly 
varying product asG at A;2 ^ e^ . The k- and /^-dependence of x^ is given by Eqs.(6,7). 

Then, going from doj^jdz to the structure function of the nucleon probed by the left- 
handed (A = L) gauge boson, Fl^x^Q"^) (see Eq.(4)), we find the soft gluon contribution to 
the P-wave component^ of 

Fl = FI + FI (17) 



denoted by F^ 



"1 dzz'^ 



Jo e^ 6t\ 
The leading contribution to F\ comes from 



mi + Q^ 



^Hereafter, the upper index in F^ corresponds to the angular monientuni L = 0, 1 of the gg'-state which 
the hght-cone W-boson transforms into. 



Therefore, Eq.(18) can be approximated by 



nM') - z:^^ ir^''-^^^Gi.„e% (20) 



2Q^ rC{^i+Q') de^asi^) ^ 
m1 + Q^Jml e^ 37r 

where C ~ 0.25. Certainly, a^^ is dominated by Q^ ^ m^. However, F^ presented in the 

form (20) allows straightforward extension to the processes induced by the charm-strange 

current, where Q^ ^ m^ and Xg ~ 2x. 

The S-wave component of F^, denoted by F]^, is very small compared to F^. Indeed, the 

^-integration in F° converges rapidly at 2; ~ 1 — ml/{mf + Q^) and 

2 2 /"^ c?zz2a5(£2) 



F°(x,Q^) ^ g^m^ / ^^G(x„5^ 
JO £ on 

G{xg,ml). (21) 



Q2 as{ml) 



Therefore, F^ < F^ and Fl ^ Fl 

The S'-wave component of the right-handed structure function 



Fn = F^^ + FA (22) 



is as follows 



JO e OTT 



'-^2 3^ G'(x„m^/2). (23) 

At Q^ < m^ the ratio (1 - zf/e'^ in Eq.(23) is flat in z and {z) ~ 1/2. 

For Q"^ ^ ml the structure function F]j is dominated by 2; ~ m'l/Q'^ and x^ ~ 2x. Hence, 

F°(a;,g^)-^^^G(a;„2mD. (24) 

The P-wave component of Fr reads 

FA(x, g^) - 2g^ /^ ^^^^^^G(x„ .^). (25) 

Jo e OTT 

For g2 < m^^ (1 - zY/e^ ^ (1 - ^)/m^ and 

FA(x,g^)-g^4^G(a:„m?). (26) 
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Figure 1: The nucleoli structure function F\ for A = S*, L, i? as a function of Q^ at a; = 5.10^^ 
in the neutrino reactions induced by the top-bottom, charm-strange and up-down quark 
currents. 



In the region of very high virtuahties of the gauge boson, CQ^ ^ m^ (see Eq.(20)), 



Jm 



~ \jr \ tXj Q ^ (Z. 



(27) 



with Xg ^ Ix. Certainly, Eq.(27) is irrelevant to the problem of a^^ : the latter is dominated by 
Q^ <^ mf. However, Eq.(27) is not entirely useless. It describes - with obvious substitutions 
- the dominant contribution to the F^ in high-Q^ processes induced by the charm-strange 
current (see Fig.l). 

From explicite expressions for F^ and Fr obtained above it is evident that for the top- 
bottom current 

Fl > Fr. (28) 

in a wide range of Q^. Another observation is that the contribution of left-handed W^-bosons 
to o"*^^ in processes induced by the top-bottom current is much smaller than that coming 
from the absorption of scalar W-bosons 



Fs^Fl 



(29) 



The scalar structure function Fs (under the name longitudinal) has been discussed in [15]. 
It was found that the higher twist corrections brought about by the non-conservation of the 
top-bottom current result in considerable enhancement of the P-wave component of 

Fs = Fl + Fl (30) 



denoted by Fg 



mf /-em de^ asie"^^ 



^^(-«^)-<TQ^il^l7"«(^»-^'- (^1) 



where e^ = mf if Q"^ < mf - ml and e^ ^ (Q^ + "^i)/4 if Q"^ > mf - ml. 

At Q^ ^ ^, -^5 is much smaller than Fg. Indeed, from Eq.(16) it follows that 



127r Jo e^ 



Fl^iL —^as{e')G{x„e'), (32) 



where Ss comes from Eq.(9) and e"^ from Eq.(ll). For Q^ ^ m^f 



2 



with S = m,h/m,t. Therefore, 

Fl{x,Q')^^^^^G{x,,m^j2) 

+ ^^Gix„ml). (34) 

Our Fs survives the limit m^f — )• oo ^. The point is that the scalar/longitudinal W-boson 
interacts with the quark current j^ = V^ — A^ which is not conserved. The vertex W ^)- tb 
which is oc d^j^ oc m^ ± m^ gives the factor m^ in the expression for Fs and cancels l/m,l 
coming from the tb quark box. Evidently, there is no real clash between Eq.(34) and the 
decoupling theorem [16]. The latter is relevant to calculating observables, for instance cross 
sections in kinematical regions well below certain heavy quark thresholds. We are dealing 
with the real top production by neutrinos of ultra-high energy, though at Q^ ^ m^. 

Two terms in Eq.(34) correspond to two very different kinematics. The first one describes 
the symmetric t6-state with uniform ^-distribution and (z) ?^ 1/2, the second one corresponds 
to the asymmetric configuration with z ^ 1 — 6'^. Both t6-states have approximately the same 
invariant mass, M^ ~ 2m^, and, consequently, probe the gluon density at nearly the same Xg. 

For Q^ ^ m,l - the substitution t -^ c suggests itself - the ratio Ss/s* is fiat in z and 
Q'^Ss/s* ~ 8 with {z) ~ 1/2 and Xg ^ 2x. Consequently, 



■^The limit ml -^ oo implies the limit of the infinite neutrino energy, 1/xg — > oo. The finite energy effects, 
built in the gluon density function , suppress the structure function in Eq. (34). Specifically, for Xg — >■ 1 the 
gluon density G is known to vanish as (1 — x^)", where n « 5 and Xg comes from Eqs.(6, 7). 
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Figure 2: The total vN charged current cross section (sohd hne) decomposed into components 
of different origin as a function of the laboratory neutrino energy. 



F°(a:,g^)^^^^g^G(x„gV4). 



(35) 



In the neutrino reactions induced by the charm-strange current the term F\ is suppressed 
at high Q^ ^ vn?^ by the factor ~ m'l/Q'^ and the scalar/longitudinal structure function 
Fs is dominated by its S-wave component, Fg [17, 18]. Specific features of the structure 
functions F\ are illustrated by Fig.l. From Fig.l and Eqs. (3,2) it follows that the left-right 
asymmetry generated by the top-bottom current is much stronger than that generated by the 
charm-strange current but the top-bottom contribution to the total neutrino nucleon cross 
section is dominated by the absorption of the scalar gauge boson. 

In our numerical studies of F\ and a'^^ we rely upon the color dipole factorization repre- 
sented by the first line of Eq.(5), where the total cross section of interaction of the th color 
dipole of the transverse size r with the nucleon target is related to the differential density of 
gluons T{xg, K^) by the equation [9] 



a{x,r) 



^ / — Z^2 "^('^ )'^(^9' *^ ^ 



K^ 



K^r"^ 



(36) 



The light cone density of the tb Fock states with the transverse size r and the fraction z of 
the W^-boson light cone momentum carried by the t-quark is [11, 12] 

2awN, r. 



I^A(^,r) 



(2vr)^ 
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(37) 



where S\ and Pa come from Eqs.(8,9). In (37) KQ^i{y) is the modified Bessel function. The 
terms proportional to KqIst) and Kf{er) describe the quark- antiquark states with the angular 
momentum L = (S-wave) and L = 1 (P-wave), respectively. The log(l/x)-evolution of 
a{x,r) is determined by the infrared regulated BFKL equation with running coupling [7, 
8]. The preferred choice of the infrared regularization gives the intercept of the pomeron 
trajectory, aip(t), in the angular momentum plane Ajp = aip(O) — 1 = 0.4 and leads to a very 
good description of the data on the proton structure functions at small x [19]. In particular, 
this means that the gluon density and all structure functions grow fast with growing 1/x, 
Fx oc x"'^"' . The analysis of the unitarity effects will be published elsewhere. The structure 
functions of the isoscalar nucleon probed by 14^-bosons of different helicity are presented in 
Fig.l. 

The charged current neutrino-nucleon cross section, o"*^^ as a function of the neutrino 
energy, E„ and its decomposition into components of different origin is shown in Fig. 2. 

Summarizing, we presented the K-factorization formulas for the differential cross section 
da\/dz(P\i which describes the absorption of scalar, left-handed and right-handed W-bosons. 
We isolated leading contributions to the related structure functions F\{x^Q'^) at low x and 
high Q^ and obtained simple and numerically accurate estimates for F\ in processes induced 
by massive quark currents. It was shown that the non-conservation of the top-bottom current 
leads to the hierarchy of structure functions Fs ^ Fl ^ Fr. Making use of the color dipole 
BFKL approach to the log(l/a;) QCD evolution we evaluated the charged current vN cross 
sections for ultra-high energy neutrino beams. 
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